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ABSTRACT 

Main-sequence stars earlier than spectral type ~ F6 or so are expected to rotate rapidly due to their radiative 
exteriors. This rapid rotation leads to an oblate stellar figure. It also induces the photosphere to be hotter (by 
up to several thousand Kel vin) at the pole th an at the equator as a result of a process called gravity darkening 
that was first predicted by Ivon Zeipel (1924). Transits of extrasolar planets across such a non-uniform, oblate 
disk yield unusual and distinctive lightcurves that can be used to determine the relative alignment of the stellar 
rotation pole and the planet orbit normal. This spin-orbit alignment can be used to constrain models of planet 
formation and evolution. Orderly planet formation and migration within a disk that is coplanar with the stellar 
equator will result in spin-orbit alignment. More violent planet-planet scattering events should yield spin-orbit 
misaligned planets. Rossiter-McLaughlin measurements of transits of lower-mass stars show that some planets 
are spin-orbit aligned, and some are not. Since Rossiter-McLaughlin measurements are difficult around rapid 
rotators, lightcurve photometry may be the best way to determine the spin-orbit alignment of planets around 
massive stars. The Kepler mission will monitor ^ 10^ of these stars within its sample. The lightcurves of any 
detected planets will allow us to probe the planet formation process around high-mass stars for the first time. 

Subject headings: techniques:photometric — eclipses — Stars:individual:Altair 
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J. INTRODUCTION 

Fifty-five transiting extrasolar planets have been discov- 
ered to date (see http://exoplanet.eu/). Most of these 
planets orbit stars that have masses near 1 Mq. The pri- 
tnary reasons for this parent star mass bias for transiting 
planets are twofold: (1) stars of spectral type later than 
K are too dim to be caught in large numbers by wide-field 
transit surveys, and (2) stars earlier than F have rotation- 
ally broadened spectral lines and inherent stellar noise that 
make high-precision radial velocity follow-up impossible at 
present. 



To address the early-star radial velocity problem, 
I Johnson et al. ( 20071) used radial velocity to survey evolved 
high-mass stars that were formerly e arly-type dwar fs when 
they were on the main se quence. Ijohnson et al. ( 2007 ) 



and Ijohnson et al 



(l2008ah described 11 known planets 
around evolved stars with 1.5 M© < M* < 3.0 Mq, 
presumably former A stars. Since then 14 new planets 
around high-mas s stars have been fou nd: NGC4349#127 b 



(Sato et al.ll2008a). HD17092b (Niedzielski et al. 2007), 14 



(|Lovis k Mavoil 12007). 81 Cetus b (|Sato et all l2008b( ). 
NGC2423#3b dLovis fc Mavod l2007l). 18 Delphinus b 



Andro medaeb (ISato et al.ll2008h), £ Aquilae b ( Sato et al . 
2008al) . HD81688b (ISato et al.l200"8ah . H D173416b (iLiu et al 
20091), HD102272b and HD102272c (iNiedzielski et al " 



2008), 6 Lyncis b (|Sato et al.ll2008bl). HD5319b (iRobinson et al 
20071), and OGLE2-TR-L9b (jSnellen et all l200^ The 



last of these, OGLE-2-TR-L9b, is the only radial velocity 
planet whose host star has A/* > 1.5 M© and lies on the 
main sequence (spectral type F3V), and this planet is the 
only one that is known to transit a massive star as well. 

The Kepler mission will discover many more transit- 
ing planets around early- type stars, if they exist. At least 
10,000 main-sequence stars earlier than spectral type F5 
should be present in the Kepler field, and will presum- 
ably be among the mission's targets. Kepler will then 
be able to characterize the distribution of planets with 
short periods around high-mass stars. These close-in plan- 
ets will complement the far-out planets recently discovered 
around A dwarfs using direct imaging (jMarois et al. I I2OO8I : 
Kalas et al.l 120081 ). and lead to a better understanding of 
how planet formation varies with stellar mass. 

Main sequence stars earlier than ^mid-F spectral type, 
those with M > 1.5 M©, are all expected to be fast rota- 
tors. The structure of these stars is such that they have 
radiative zones in their outermost layers, instead of a con- 
vective zone near the surface like for our Sun. The ex- 
terior convection in later-type stars drives surface mag- 
netic activity, which in turn drives strong stellar winds that 
sap the star's angular momentum with time. Early-type, 
exterior-radiative stars retain their youthful high angular 
mo menta, with some sp inning at near the breakup speed 
(e.g. [Hansen et ah 2004 . and references therein). As a re- 
sult, early- type main sequence stars can be significantly 
oblate. 
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The rotation induces an equator-to-pole gradient in 
the effective acce leration due to gravity g at the surface. 
von Zeipe showed that in such a case the tempera- 

ture of the star varies from equator to pole as well, a phe- 
nomenon called gravity-darkening. The Von Zeipel Theo- 
rem thus predicts that the flux emitted from the surface 
of a rapidly rotating star is proportional to the local effec- 
tive gravity. Thus the effect induces cooler temperatures 
(and hence lower emitted fluxes) at a star's equator, and 
hotter temperatures at the poles. The basic predictions 
of von Zeipel theory were dramatically confirmed by re- 
cent optical interferom etric observations of Vega (a Lyrae) 
( Peterson et al.ll2006b | . explaining residua ls in earlier near- 
IR interfero metry by fciardi et al.l ( 200lL ) and Altair (a 
Aquilac) (M onnier et al.l 20071 ). Gravitv darkening is used 
regularly to character ize close binary star s ystems from 
their lightcurves {i.e. iDiurasevic et al. I l2003h . In binary 
systems that interact gravitationally, tides can also reduce 
the effective gravity, resulting in gravity darkening. 

If planets around fast-rotating stars formed in-situ from 
the protostellar disk or migrated to their present locations 
within that disk, then those planets might be expected to 
orbit near their stars' equatorial planes. If those planets 
transit, then their orbital inclinations i are near 90° (using 
radial velocity teams' definition of i as the angle between 
the planet's orbit pole and the plane of the sky). The stel- 
lar orbit pole would then be nearly coincident with the 
planet's orbit pole, giving a stellar obliquity tp of near 0°. 
In transit, such a planet's chord across its star's disk would 
be perpendicular to the projected stellar rotation pole. In 
this paper, I call these "spin-orbit aligned" planets. If 
planets have experienced planet-planet scattering events in 
their past, however, they might be expected to show sig - 
nificant spin-orbit misalignment ( Juric fc Tremaindl2008l ). 



Rossiter-McLaughlin (R-M) measurements of the ra- 
dial velocity of slowly-rotating stars during planetary 
transits have been highly successful at determining spin- 
orbit alignments. Relatively low-precision R- M measure- 



ments of s pin-orbit alignement w ere made by iPont et al . 
( 20091) andlMoutou et all (|2009l ). IWinn et alJ ([20061) and 



Winn et al.((|2007t ) made early measurements of HD189733 
and HP 14750 6 , show ing them to be spin-orbit aligned. 
Johnson et al. I (|2008bl) determined th a t HAT -P-1 is nearly 
aligned (3.7° ± 2.1°). iNarita et all (|2009l ) showed that 
HD17156 is nearly spi n-orbit alig n ed, w ith a misalign- 



ment of 10° ± 5.1° 



Wolf et al. 



(|2007n determined a 
sinnlarly misalignment of 12° ± 15° in the HD149026 
system, con sistent with sp i n-orb it alignment. On the 
other hand, iHebrard et all (|2009t ) show a striking spin- 
orbit misalignment of 70° ± 15° in the X O-3 pl anetary 
system . Highly precise measurements from iTriaud et al.l 
(|2009f ) show a tiny but significant spin-orbit misalignment 
of 0.85° ± 0.3° in the CoRoT-3 system. Unfortunately, 
Rossiter-McLaughlin measurements will likely be much 
more challenging rapidly-rotating stars due to their high 
inherent radial velocity noise. 

Because more massive stars rotate much faster than the 
Sun, the transit lightcurves for the planets that Kepler 



will discover around them will be qualitatively and quan- 
titatively differe nt from those for planets orbiting slowly- 
rotating stars. Russelll ( 1939f l considered this effect for 
eclipsing binary stars, and here I consider the effects for 
transiting planets. Stellar oblateness will alter the times of 
transit ingress and egress, and the overall transit duration, 
somewhat complementary to the effects of o blate planets 
(jBarnes fc FortnevI 120031 : ISeager fc Huill2002h . When the 
star's spin pole and the planet's orbit pole are ahgned, 
the Von Zeipel effect will cause systematic errors in ra- 
dius determinations for the star and the planet, and will 
lead to broadband color variations during transit. If the 
stellar spin pole and planetary orbit pole are not aligned, 
then bizarre transit lightcurves result that can be used to 
constrain both the stellar spin pole direction and the spin- 
orbit alignment. In this paper I investigate the effect that a 
fast-rotating star has on the lightcurves of transiting extra- 
solar planets in preparation for the results expected from 
Kepler. 

2. SYNTHETIC TRANSIT LIGHTCURVES 
2.1. Algorithm 

In order to generate synthetic transit lightcurves with 
fast-rotating stars, I modified t he alg orithm originally 
developed for Barnes fc FortnevI (2003) and extended in 
Barnes fc FortnevI (|2004l ) and lBarn'e^(|2007l l. The algo- 
rithm numerically integrates the total flux coming from 
the uneclipsed star, Fq, in polar coordinates centered on 
the projected center of the star in the plane of the sky such 
that: 



Fa 



I{r, 0) de dr 



(1) 



where Rcq is the radius of the star at its equator (see Fig- 
ure [T] for a schematic of some of the geometric variables), 
r and 9 are measured from the stellar center and counter- 
clockwise from the x-axis respectively, and /(r, 6) is the 
star's intensity at point {r,9). It then evaluates the appar- 
ent stellar flux at time t, F{t), relative to the out-of-transit 
flux Fq, by subtracting the amount of stellar flux blocked 
by the planet from _Fo- 



blocked 



T{r,0,t)I{r,6) d9 dr , (2) 



and 



blocked 



Fo 



(3) 



where r{r,6,t) = 1 if the planet is blocking starlight at 
position r, 9 and time t, and r(r, 9,t) = if not. 

The difference in the case of fast-rotating stars is that 
the rotation induces these stars' equators to bulge out- 
ward. Hence the polar integral in Equations [T] and [2] no 
longer properly accounts for the symmetry of the prob- 
lem. I quantify this effect here as the star's oblateness, /, 
defined to be / = ^"''jj.^^"'" where i?poio is the star's ra- 
dius along the rotational pole. I assume throughout that 
the star's resulting shape can be considered a MacLaurin 
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spheroid. The value relevant for the integrals in Equations 
[1] and O though, is the effective oblateness /oft, which I 
define to be the apparent oblateness of the star when pro- 
jected into the plane of the sky. The effective oblateness is 
related to the actual oblateness by the stellar obliquity 
where ip — if the stellar rotation axis resides in the plane 
of the sky: 



/off = 1 - •^(1-/)^ 



is^ (f + sin^ Lp 



(4) 



The relationship is not a simple cosine owing to three- 
dimensional geometry. I will show how this expression can 
be derived a few paragraphs below. 

With /eff in hand, it becomes straightforward to incor- 
porate the stellar asphericity. In order to avoid complex 
and computationally-intensive elliptical integrals, I execute 
a substitution for r and 6 in Equations [1] and [21 I instead 
choose to integrate over r' and 9' , where r' and 9' are cho- 
sen so as to 'pop' the star into a spherical shape in r'-d' 
space. To do this, I first convert the true projected r and 
6 measured from the star's center into x = r cos 9 and 
y = rsin^. I assume that the projected stellar rotation 
axis is parallel to the y-axis for simplicity - the true orien- 
tation will not be known, in general, but does not matter 
since it does not affect the measured stellar flux. I then let 
x' = X and 



y 



y 



(1 - /off) 



(5) 



and then set r' = y/ x'"^ + y'^ and 9' — atan2(i;', x') where 
atan2 is the C computer language arctangent function that 
returns a true 4-quadrant-capable angle from x and y val- 
ues. The substitution in Equation [5] works equally well if 
you were to choose to integrate the stellar flux in cartesian 
xy coordinates instead of the polar integral that I use. 



Equations [T] and [5] now become: 



i^o = (1 - /cff) 



2tt 



Iir',9') d9' dr' 



(6) 



and 



Fblocked{t) — (1 — ./off) 



2Tr 



V{r' ,9' ,t)I{r' ,9') d9' dr' . 

(7) 

The (1 — /off) factor is introduced by the coordinate trans- 
formation: 'popping' the oblate star out into a circle over- 
estimates its projected area, and hence the emitted flux, by 
(1 — /eff)~^- In the end the factor is irrelevant. When the 
results are plugged into Equation [31 it drops out. Hence 
the lightcurve generation algorithm as implemented does 
not use the factor (1 — f^s) explicitly at all. 

All that's left then is to determine I{r',9') and then 
to integrate it. Obtaining /(r, 9) is straightforward but 
nontrivial. I first break out the stellar limb darkening from 
the normal emission, and assume blackbody radiation: 



Iir',9') = Bx{T{r',9'))Lir',9') 



(8) 



where B\(T) is the blackbody function (or your desired 
stellar emission as a function of temperature at a given 



wavelength), T{r',9') is the temperature at a given point 
on the stellar disk, and L{r',9') is the stellar limb dark- 
ening at that point. It would be possible to incorporate a 
more realistic stellar emission spectrum rather than to as- 
sume it to be a blackbody, but since rapidly-rotating stars 
are mostly of early spectral types, the differences are not 
significant at the level of the present investigation. 



A modified version of the Von Zeipel theorem (jvon Zeipel 



19241 ) determines the stellar temperature T at every point. 
The temper ature on . the s urface of a rapidly rotating star 
is given by (|Maedeill2009l ) 



T = 



9' 

pole p 
fl'pole 



(9) 



where g is the magnitude of the local effective surface grav- 
ity and Tpoic and gpoic are the pole's temperature and sur- 
face gravity, respectively. The value /3 is known as the 
gravity-darkening parameter. Its nominal value is 0.25 for 
a purely radiative star, as derived by von Zeipel; for our 
numerical calculations I use an empirically measured /?, as 
detailed below. The local surface gravity vector g has two 
terms; one Newtonian, and one centrifugal: 

CM 

9^-^f + ^'R±rl (10) 

where G is the universal gravitational constant, M, is the 
stellar mass, J7 is the stellar rotation rate in radians per 
second, R is the distance from the star center to the point 
in question, and R± is the distance from the star's rota- 
tion axis to the point in question. The two hatted sym- 
bols, f and rj^ are unit vectors pointing to the point in 
question from the stellar center and stellar rotation axis 
respectively. 

So, then, to know / you need to know g, and to know g 
you need to know the three-dimensional vector position of 
each point that you see on the star. In this case you already 
know r' and 9' , from which you can trivially derive x and 
y, the location of 2-d projection in the plane of the sky of 
the point of interest with respect to the center of the star. 
What is left then is to determine z. This is nontrivial. 

The geometrical constraint on z is that its value must 
conform to the surface of an oblate spheroid given x, y, 
the stellar radius i?eq, and the oblateness /. This is easy 
enough when the stellar obliquity is zero, and thus when 
the y-axis is parallel to the stellar rotation axis. So I define 
a new set of coordinate axes with the same origin as the 
x-y-z system, at the center of the star. However this new 
system is rotated in the y-z plane (i.e. around the x-axis) 
by an angle (p, the star's obliquity to the plane of the sky. 
Call this the xq, yo, zq system, where 



Xq 

yo 



y cos{(p) + z sin{(p) , and 
—y sm{ip) + z cos{ip) . (11) 



In this new obliquity-rotated coordinate system, the sur- 
face of the star's photosphere follows 



yi 



+ zl 



Rt 



(12) 
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Plugging the definitions of xq, yo, and zq from Equations 
111! into Equation 1121 I solve for z in terms of a known x 
and y. The solution to the resulting quadratic is 



'2y{l — (1 — /)^) simp cos ip + \fd 



(13) 



(14) 



2((1 - /)2cosV + sin^ </') 
where the determinant d is 

d= V(i _ (1 - /2))2sin2(^cosV + 
-4((cos2(p(l - /)2 +sin2 ^) X 
((y2 sin^ ^ - R\ + x2)(l - f) + y2 cos^ ^)) . 

I choose the positive root of the determinant as the neg- 
ative root represents the invisible second interception of 
the line-of-sight with the photosphere that occurs on the 
back side of the star as seen from Earth. Equation |4] above 
can be derived by setting the determinant d equal to zero 
to establish the outer edge of the star's disk as seen from 
Earth, and then solving for the proper /eff to reproduce 
that disk. 

I now have all of the necessary parameters to compute 
the flux coming from each point on the star. To do so, at 
each X, y point use Equation 1131 to get z, and then plug 
z into Equations [TT] to get xq, y^, and zq. Do the vector 
addition to get g from Equation [10] in xq, j/Oj space 
where 



R 



-^IQ + -^fco , 

rL± -fix 



R = 



R, = 



Xq + Vo + ^0 



Then plug g = \g\ into Von Zeipel's Equation (EquationlHl) 
to get T, and then derive a flux from T using a blackbody 
curve or your choice of a more sophisticated emitted flux. 

2.2. Parameters 

In order to generate appropriate and representative 
lightcurves that can be used for comparison to transits 
yet undiscovered, I calculate all transit lightcurves as if 
the parent star were Altair (a Aquilae). The true host 
stars for Kepler-detected transiting planets will show vary- 
ing stellar masses, polar temperatures, radii, and rotation 
rates. I elect to use Altair because I think that its spec- 
tral type (A7V) is broadly representative of the majority 
of the expected fast-rotating stars in the Kepler sample, 
and because its parameters are well-chara cterized by inter- 
ferometric imaging (Monnier et al.l 120071 ). Not all sets of 
parameters produce physically plausible stars; I avoid non- 
physical combinations by only using this one set of known 
stellar values. 



The s pecific stellar parame ters that I use are M, 
1.8 Mq ((Peterson et al.ll2006al ). R^q 



■0 1 ^polc 



j |2006a[ ). i?eq = 2.029 Rg 
8450 K, /3 = 0.190, / = 0.1947, and a stellar rotation pe- 
riod of 8.64 hours (all as directly measured for Altair by 



Monnier et al.ll2007l ). For the planet, I assume Rp = Rjup 
and an orbit semimajor axis of 0.05 AU (corresponding to 
a period of 3.04 days) for familiarity with the lightcurves of 
known transiting hot Jupiters. Transit lightcurve shapes 
are invariant with orbit period; only the timescale changes. 
Hence the curves that I show here can be converted for 
different-period planets by stretching the x-axis. 

In fitting I adjust for the Ci an d Ca lim b darkening co- 
efhcients outlined in iBrown et al.l (2001). I generate the 
synthetic hghtcurves using ci = 0.640 and C2 = 0.0. I 
also assume a monochromatic observation at 0.51 microns 
wavelength except where otherwise noted. 

3. SPIN-ORBIT ALIGNED 

I show synthetic transit lightcurves of planets in spin- 
orbit aligned geometries in Figure [51 for various transit 
impact parameters b. The lightcurves are symmetrical. 
Furthermore, since the parts of the star that the planet 
transit chord passes over all have the same temperature 
because they equidistant from the stellar rotation pole, 
the transit bottom shows normal limb-darkening curva- 
ture. As a result, the specific lightcurve shapes for spin- 
orbit aligned transiting planets are indistinguishable from 
those for planets orbiting spherical, slow-rotating stars. 

The transit depths and durations are different, though. 
In particular, the transit depth increases with impact pa- 
rameter, opposite of the case for spherical stars. If you 
were to fit these lightcurves assuming a spherical star, then 
the resulting best-fit transit parameters would be differ- 
ent from the actual ones. In order to show this effect, I 
fit the synthetic transit lightcurves show in Figure [2| using 
a Le avenberg-Marquard t chi-squared mi nimization scheme 
from [Press et al.l ( 1992[ ) as described in Barnes fc Fortnev 
(|2003ir and show the resulting best-fit parameters in Ta- 
ble [S 

For planets transiting across the center of the star 
{b = 0.0), the best-fit parameters retain easily-interpreted 
astrophysical meaning. The ratio of radii Rp/ R^, comes 
from the total transit depth, which in this case is nearly 
identical to that for the case of a 6 = 0.0 transit across 
a spherical star with radius equal to the equatorial radius 
of the fast-rotating star. In that case, while using the Al- 
tair parameters that I use here, the total larger projected 
area of the stellar disk makes up for lower flux coming 
from the polar regions, yielding similar net stellar fiux. 
The fitting algorithm thus gets Rp/R^, correct, assuming 
i?* = Req and not the real average radius of the projected 
disk, (1 — ^)Req. The total duration of the transit and the 
duration of ingress and egress fix the impact parameter, 
which the fit correctly determines to be near 0.0 {i ^ 90°), 
and the stellar radius, which is very close to the star's true 
equatorial radius. Hence the planetary radius also comes 
out correctly. Limb darkening is right because of the same- 
temperature effect described above. 

As I consider higher and higher impact parameters, 
shown in Figure [2] as fractions of Altair's polar radius 
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Table 1 

Best-fit transit parameters assuming spherical star for various transit geometries. 



Stellar 


Planet's 


Transit 


Best- 


fit 


Parameters 






Obliquity 


Azimuth 


Impact 








Limb 


Limb 


(V', °) 


(«, °) 


Parameter 


R* 


Rp 


Inclination 


Dark 1 


Dark 2 






(6, Rpole) 


R© 


Rjup 


ii, °) 


Cl 


C2 




0° 


0° 


0.0 Rpoic 


2.050 R© 


0.999 Rjup 


88.23° 


0.612 




0° 


0° 


0.3 Rpoic 


2.014 Rq 


1.034 Rjup 


86.95° 


0.621 




0° 


0° 


0.6 Rpoio 


1.970 Ro 


1.141 Rjup 


84.07° 


0.626 




0° 


0° 


0.9 Rpoic 


1.937 R© 


1.291 Rjup 


80.94° 


0.582 




30° 


0° 


0.9 Rpoic 


1.91 Rq 


1.27 Rjup 


81.70° 


0.826 


1.000 * 


30° 


0° 


—0.6 Rpoic 


2.02 Rq 


1.24 Rj„p 


84.01° 


0.834 


1.000 * 


30° 


0° 


0.3 Rpoic 


2.08 Rq 


1.16 Rjup 


85.93° 


0.805 


1.000 * 


30° 


0° 


0.0 Rpoic 


2.04 Rq 


1.02 Rj„p 


88.70° 


0.723 


1.000 


30° 


0° 


0.3 Rpoic 


2.02 Rq 


1.01 Rjup 


89.97° 


0.726 


0.859 


30° 


0° 


0.6 Rpoic 


2.123 Rq 


0.987 Rjup 


83.05° 


0.465 




30° 


0° 


0.9 Rpoic 


2.143 Rq 


1.059 Rjup 


80.19° 


0.268 




90° 


0° 


0.0 Rpoic 


2.30 Rq 


1.17 Rjup 


80.30° 


0.916 


1.000 * 


90° 


0° 


0.3 Rpoic 


2.06 Rq 


1.00 Rjup 


84.63° 


1.000 


1.000 * 


90° 


0° 


0.6 Rpoic 


2.05 Rq 


1.01 Rjup 


87.10° 


1.000 


1.000 * 


90° 


0° 


0.9 Rpoic 


2.04 Rq 


1.01 Rjup 


89.47° 


1.000 


1.000 * 


0° 


90° 


0.0 Rpoic 


1.89 Rq 


1.12 Rjup 


84.57° 


-1.000 


-1.000 * 


0° 


90° 


0.3 Rpoic 


1.91 Rq 


1.11 Rjup 


84.00° 


-1.000 


-1.000 * 


0° 


90° 


0.6 Rpole 


1.92 Rq 


1.04 Rj^p 


82.94° 


-1.000 


-1.000 * 


0° 


90° 


0.9 Rpole 


1.85 Rq 


0.91 Rjup 


82.12° 


-0.585 


-1.000 * 


Note. — 


Tabled here 


are the best-fit parameters i?,, Rp, i, ci 


and C2 for fits 


assuming 


; spherical 



stars of synthetic lightcurves of hypothetical Altair-planet systems with Rp = Rjup and the planet 
at 0.05 AU. The stellar obliquity to the plane of the sky ip for each synthetic curve is listed at left, 
followed by the angle a between the planet's orbit pole and the stellar rotation axis projected into 
the plane of the sky. When fixing C2 at zero and fitting only for ci resulted in a high-quality fit to 
the data, C2 is not listed. Fits that have a '*' next to their value for C2 did not produce good fits even 
with the second limb darkening parameter. 
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Time From Inferior Conjunction (s) 



Fig. 2. — Synthetic lightcurves for transiting 1 Rjup in a spin-orbit aligned 0.05 AU orbit around an Altair-like star are 
plotted. The four curves correspond to planets with transit impact parameters of 6 = 0.0 Rpoie (solid), h ~ 0.3 Rpoie (dashed), 
b = 0.6 Rpoie (dot-dashed), and b = 0.9 Rpoie (dotted). The curves' shapes are indistinguishable from transits of slow-rotating, 
spherical stars using different parameters. Transits nearer to the pole are deeper because the stellar photosphere is hotter 
there due to the Van Zeipel effect. 



1.63 Rq, the best-fit parameters deviate further and fur- 
ther from the actual values. As the planet transits across 
hotter, brighter parts of the stellar photosphere near the 
pole, the best-fit value for planet radius increases by up 
to 30% above the input planet radius. The stellar radius, 
driven by the somewhat longer transit durations, drops 
slightly but stays closer to the input star's equatorial ra- 
dius than it does to the true average radius of the stellar 
disk. Less importantly, the inclination and limb darkening 
parameter are somewhat underestimated. 

The end result of these calculations is that if the planets 
orbiting fast-rotating stars are spin-orbit aligned, it may 
not be evident at first glance. The transit lightcurves will 
not stick out. Instead the measured transit parameters' 
errors will be systematic in nature, and the fits will still be 
good. 

If the spin-orbit aligned planet case turns out to be 
prevalent around fast-rotating stars in the Kepler sample, 
then measuring the spin-orbit alignment for planets will 
require a separate identification of the star's fast-rotating 
status. Kepler'' s photometric precision should be good 
enough to identify the star's rotation rate in long time- 
series photometry. If the planets transiting those stars 
show normal-looking lightcurves, that would indicate that 
the planet is spin-orbit aligned. If most or all planets are 
spin-orbit aligned, then they likely either formed in-situ or 
migrated there, and were not scattered by close encounters 
with other planets. 

4. SPIN-ORBIT MISALIGNED 

If instead the transiting planets around fast-rotating 
stars experienced close encounters that threw them out 



of spin- aligned orbits, then more exotic lightcurves result. 
Recent Rossiter-McLaughlin results from transiting plan- 
etary systems indicate that this si tuation may be mor e 



common than previously thought ( Hebrard et al. 20091 ) . 



Given the difficulty of radial velocity, and thus Rossiter- 
McLaughlin, measurements around early- type stars, tran- 
sit photometry may then be the best way to measure the 
spin-orbit alignment in those systems. 

4.1. SYMMETRIC 

Just because the transit lightcurve of a planet around 
a fast-rotating star is symmetric does not mean that the 
planet is spin-orbit aligned. Any planet whose transit 
chord is perpendicular to the projected stellar rotation 
pole, i.e. for which the angle between the transit chord 
and the projected pole, a, is 0.0°, will show a symmetric 
lightcurve. I plot the lightcurve for a few such hypothetical 
planets in Figures [3] and ID Figure [3] shows planets around 
a star with obliquity </? = 30°. Figure |4] shows transits for 
a star viewed pole-on with if — 90°. 

Although these lightcurves are symmetrical about the 
mid-transit point, unlike the spin-aligned case they are 
not well-fit using a model that assumes a spherical star. 
In the if = 30° case the stellar obliquity removes the sym- 
metry in the stellar disk, thereby making transits toward 
the north and south poles distinct from one another. I ar- 
bitrarily choose to define those transits toward the north 
of the center of the stellar disk to have a negative impact 
parameter so as to differentiate the two possible cases. 

In the positive impact parameter cases in Figure [3l the 
planet crosses the cooler equatorial regions of the star. 
Hence those transits are relatively shallow. The geom- 
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Fig. 3. — Synthetic lightcurves for transiting 1 Rjup in a 0.05 AU orbit around an Altair-like star with obhquity 30° are pfotted. 
The seven curves correspond to planets with transit impact parameters of 6 = —0.9 Rpoie (thick dotted), b = —0.6 Rpoie (thick 
dot-dashed), b = —0.3 Rpoie (thick dashed), b = 0.0 Rpoie (thick sohd), b = 0.3 Rpoie (dashed), b = 0.6 Rpoie (dot-dashed), 
and b — 0.9 Rpoie (dotted). The curves' shapes are start to differ from those of transits of slow-rotating, spherical stars. The 
ingress and egress of transits at opposite impact parameter are nearly the same. The transit bottoms differ. Transits near 
the hot stellar north pole are deep and show severe curvature of the transit bottom. 




Time From Inferior Conjunction (s) 

Fig. 4. — Synthetic lightcurves for transiting 1 Rjup in a 0.05 AU orbit around an Altair-like star with obliquity 90° (pole-on) 
are plotted. The four curves correspond to planets with transit impact parameters oi b = 0.0 Rpoie (solid), b — 0.3 Rpoie 
(dashed), b = 0.6 Rpoie (dot-dashed), and b = 0.9 Rpoie (dotted). The lightcurves show a pronounced 'U'-shape that is 
typically considered to be characteristic of grazing eclipsing binary stars rather than planets. The extreme curvature of the 
usually flat planet transit bottom makes the locations of second and third contact difficult to discern. Due to symmetry 
around the stellar pole, all transits with these same impact parameters will look the same regardless of a. 
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Fig. 5. — Synthetic lightcurves for transiting 1 Rjup in a 0.05 AU orbit around an Altair-like star with obhquity 0° (equator- 
on) are plotted, similar to Figure [2l but this time with an azimuthal angle of a = 90° . The four curves correspond to planets 
with transit impact parameters of 6 = 0.0 Rpoie (solid), b = 0.3 Rpoie (dashed), b = 0.6 Rpoie (dot-dashed), and b = 0.9 Rpoie 
(dotted). This unlikely 90° transit geometry also produces symmetric transit lightcurves, albeit highly unusual ones. These 
curves are deepest near second and third contact, and shallow at mid-transit. If assuming a spherical-star model, then the 
interpretation of this "double-horned" structure might be a negative limb-darkening coefficient possibly associated with a 
temperature inversion in the stellar atmosphere; however, such a model fits these data very poorly (see text). 



etry of the orthographic projection of the star as seen 
from Earth means that for the high impact parameter, 
b = 0.6 Rpoie and b = 0.9 Rpoie cases, the planet cov- 
ers more southerly stellar latitudes at ingress and egress 
than it does at mid-transit. Because those more southerly 
latitudes are hotter, the net result is a flatter lightcurve 
bottom than occurs for spherical stars with the same limb 
darkening parameter (ci = 0.640). The flatness appears as 
a lower limb darkening parameter ci in Table [T] For the 
b = 0.6 Rpoie and b — 0.9 Rpoie cases the model overesti- 
mates the star's radius; the best-fit radius is greater than 
even the star's equatorial radius. 

For b — 0.0 Rpoie and b — 0.3 Rpoie in the (p — 30° 
case, the planet transits closer to the star's north pole, 
and hence the transit depth becomes progressively greater. 
The ingress and egress are again at more southerly lati- 
tudes than mid-transit. But now mid-transit is north of 
the equator, so the ingress and egress are cooler than mid- 
transit, enhancing the curvature of the lightcurve bottom. 
I was unable to find satisfactory fits t o either of the s e syn - 
thetic lightcurves using just the first iBrown et al.l (|200ir ) 
limb darkening coefficient ci alone, even though it was the 
only one used to generate the synthetic lightcurve. How- 
ever, by allowing the fitting algorithm to optimize C2 as 
well I managed to find reasonable fits for b = 0.0 Rpoie and 
b = 0.3 Rpoie- 

When fitting b = —0.3 Rpoie, b = —0.6 Rpoie, and 
b = —0.9 Rpoie not even two-limb-darkening fits were ac- 
ceptable. The curvature on these transits is so extreme 
that even though they are symmetric, they are discern- 



able from reasonable spherical-star models because of the 
severity of their curvature. In fact, the curvature of these 
planets' transit bottoms are so severe that the lightcurves 
look 'V'-shaped instead of 'U'-shaped. The transit bottom 
is difficult to discern. The resulting rounded shape might 
be mistaken for an eclipsing binary star. 

Fits for transits of pole-on oriented fast-rotating stars 
are similar. As shown in Figured! the center of the stellar 
disk is particularly bright relative to the limb owing to the 
center being the hot pole and the limb being the cool equa- 
tor. Thus here again the curvature of the transit lightcurve 
between second and third contacts is large, and cannot be 
satisfactorily fit by a spherical-star model. In both these 
cases the difficulty in identifying these type of planetary 
transits in the Kepler data will probably be recognizing 
that they are not eclipsing binaries, not in thinking that 
they are planets around spherical stars. 

An observer would also measure a symmetrical transit 
lightcurve if a planet orbits a zero obliquity {ip — 0.0°) 
star such that its orbit plane contains the stellar rotation 
pole, i.e. with a — 90°. I show synthetic lightcurves for 
such a situation in Figure [S] These are strange transits. 
By moving parallel to the star's orbit pole, the planets first 
encounter the limb-darkened edge of the star. Shortly after 
second contact, they cover the brightest part of the stellar 
photosphere underneath the transit chord. The lightcurve 
is deepest there. Over the equator at mid-transit the tran- 
sit depth is shallow, and then the process repeats itself 
backwards on egress. 
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The resulting lightcurve has a "double-horned" struc- 
ture that bears a resemblance to what a transit across a 
limb-brightened star would look like. However, spherical- 
star model fits cannot reproduce the specific structure. In 
particular, a limb-brightened star would have sharp points 
at second and third contacts where the derivitive is discon- 
tinuous and changes sign. The fast-rotating star ip — 0° 
a = 90° transits in Figure [5] instead have rounded horns at 
second and third contact owing to the combination of limb- 
and gravity-darkening. The resulting rounded double-horn 
shape is both characteristic and diagnostic of this kind of 
transit. 

4.2. ASYMMETRIC 

If the stellar orientation and transit chord have more ex- 
otic geometries, then highly unusual asymmetric lightcurves 
result. A star with obliquity ip — 30° and transit chord 
azimuth a = 90° (motion parallel to the projected stellar 
rotation axis) presents an easy such case to visualize (Fig- 
ure O. Because the stellar disk presents only a left-right 
symmetry, these transits pass over very different photo- 
spheric temperatures in the first half of the transit relative 
to the second half. 

In the cp = 30°, a — 90°, impact parameter b = 0.0 Rpoic 
case, the planet first crosses a darkened stellar limb be- 
fore occulting the star's hot north pole. The resulting 
lightcurve is deepest at this point. After mid-transit the 
planet occults the cool equator for a shallower transit 
depth, but then the depth increases again near third con- 
tact as the planet covers higher southern latitudes near the 
southern stellar limb. The same process occurs in a more 
abbreviated fashion as the impact parameter b increases. 

If the transit chord azimuth is oblique, then the 
lightcurves can become quite complex as shown in Fig- 
ure [71 The central 6 = 0.0 Rpoic transit resembles those 
from Figure [6l The more northerly transits (those defined 
to have negative impact parameters) also show strong 
lightcurve asymmetries, with deep first halves and shal- 
lowing second halves. Some of the lightcurves turn over 
again near third contact, but the ones with more negative 
impact parameters do not. 

For postitive, more southerly impact parameters, the 
depth asymmetry decreases. Futher from the hot pole 
the temperatures under the transit chord are more uni- 
form. The b = 0.9 Rpoic transit is particularly interesting. 
It shows a nearly uniform depth in time, but while the 
ingress is long and gently curving, the egress is abrupt. A 
higher photometric precision would be required to defini- 
tively identify such a transit lightcurve as being one from 
a fast-rotating star than would be necessary for some of 
the more spectacularly asymmetric lightcurves. 

With oblique azimuths also come uncentered lightcurves. 
Because of the stellar oblateness, the center of these tran- 
sits does not occur at the time of the planet's inferior con- 
junction, that is, when the planet is closest to the Earth. 
The total discrepancy can be up to a few tenths of the 



total transit duration, depending on the stellar oblateness 
and the transit geometry. 

The oblique azimuth transits of fast-rotating stars also 
introduce an asymmetry in the duration of transit ingress 
and egress. This effect comes about because the angle 
between the stellar limb and the transit chord is different 
on ingress and egress, thereby causing the transiting planet 
to take different amounts of time to cross the limb in each 
case. 

Lightcurve asymmetries can occur in zero-obliquity 
{ip — 0°) stars as well when the transit chord azimuth 
is oblique. I show such a case in Figure [H with ip = 0° and 
a = 60°. The central b = 0.0 Rpoic transit is symmetric, 
similar to those from Figure [5l Non-central transits are 
asymmetric and resemble those from Figure [6] with one 
half of the transit being deeper than the other. In this 
ip = 0° oblique case, the transits with negative impact pa- 
rameter arc the time-inverse of those with positive impact 
parameter. 

5. COLOR EFFECTS 

Due to the temperature nonuniformity across the stel- 
lar disk, transit lightcurves around rapidly-rotating stars 
differ as a function of wavelength. I show an example of 
this effect in Figure [9l While a hotter blackbody radiates 
more flux at all wavelengths than a colder one, the flux 
ratio (hotter over colder) is greater at the Wien end of the 
blackbody curve. Hence the effects of rapid rotation in 
transit lightcurves are maximized at short wavelengths. 

At wavelengths much longer than the blackbody emis- 
sion peak, the flux ratio between the hotter and colder 
areas approaches the temperature ratio in Kelvin. Uti- 
lizing multiple wavelength lightcurve observations would 
strongly constrain the temperature structure on the stel- 
lar disk in addition to the transit parameters. This effect 
will be important when fltting for a fully-consistent model 
for the star (instead of assuming Altair's parameters as I 
do in this paper). 

6. CONCLUSION 

The gravity-darkeni ng effect for rapidly rotating, oblate 
stars first predicted bv lvon Zeipell (|l924) allows for highly 
unusual lightcurves when such a star is transited by an 
extrasolar planet. The distinctive lightcurves allow for 
a solely photometric determination of the relative align- 
ment between the stellar spin axis and the planet's orbit 
normal (spin-orbit alignment), a measurement that usu- 
ally requires radial velocity measurements of the Rossiter- 
McLaughlin effect. The alignment bears a fingerprint of 
the planet's formation and evolution. Planets that form in 
an orderly fashion and migrate within a disk ought to end 
up coplanar with the stellar equator. Those that have ex- 
perienced planet-planet scattering events ought in general 
to not be spin-orbit aligned. 

Spin-orbit aligned planets around fast-rotating stars 
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Fig. 6. — Synthetic lightcurves for transiting 1 Rjup in a 0.05 AU orbit around an Altair-like star with obhquity 30° are 
plotted, similar to Figure [3l but this time with an aziniuthal angle of a = 90°. The four curves correspond to planets with 
transit impact parameters of 5 = 0.0 Rpoie (solid), b — 0.3 Rpoie (dashed), b — 0.6 Rpoie (dot-dashed), and b — 0.9 Rpoie 
(dotted). The resulting lightcurves are highly asymmetric, being deeper on the side of the lightcurve where the planet passes 
over the hot northern stellar pole. 




Fig. 7. — Synthetic lightcurves for transiting 1 Rjup in a 0.05 AU orbit around an Altair-like star with obliquity 30° are 
plotted, similar to Figure [3] but an azimuthal angle of a = 60° . The seven curves correspond to planets with transit impact 
parameters of 6 = —0.9 Rpoie (thick dotted), b = —0.6 Rpoie (thick dot-dashed), b = —0.3 Rpoie (thick dashed), b — 0.0 Rpoie 
(thick sohd), b = 0.3 Rpde (dashed), b = 0.6 Rpoie (dot-dashed), and b — 0.9 Rpoie (dotted). With this oblique azimuth, the 
photometric lightcurve center does not correspond with the point where the planet is nearest to Earth. The lightcurves show 
a diversity of complex asymmetric shapes as a function of impact parameter. 
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Fig. 8. — Synthetic lightcurves for transiting 1 Rjup in a 0.05 AU orbit around an Altair-like star are plotted, similar to Figure 
[7] but with an obliquity (p — 0° . The seven curves correspond to planets with transit impact parameters of 6 = —0.9 Rpoie 
(thick dotted), b = —0.6 Rpoie (thick dot-dashed), b = —0.3 Rpoie (thick dashed), b = 0.0 Rpoie (thick soUd), b = 0.3 Rpoie 
(dashed), b = 0.6 Rpoie (dot-dashed), and b — 0.9 Rpoie (dotted). Due to the symmetry of the stellar disk in this case 
the b = 0.0 Rpoie case shows a symmetric lightcurve, and positive impact parameters have time-reversed lightcurves when 
compared to their negative counterparts. 
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Fig. 9. — Synthetic lightcurves for transiting 1 Rjup in a 0.05 AU orbit around an Altair-like star with (p = 60° are plotted. 
The 4 curves all correspond to planets with a transit impact parameter of 6 = —0.3 Rpoie- The different curves correspond 
to what a lightcurve would look like if it were acquired at different wavelenghts: 0.25 ^m (solid line), 0.375 /xm (dotted line). 
Note that wavelength is all that is changed here - in particular I have employed the identical limb darkening parameter 
of ci = 0.64 for each wavelength. Though this is unphysical in that the limb darkening should diminish with increasing 
wavelength, this way makes it easiest to identify just the differences due to wavelength changes. The contrast on the stellar 
disk is greater on the Wien (short wavelength) end of the blackbody curve, while the contrast is muted on the Rayleigh-Jeans 
(long wavelength) side. 
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show lightcurvcs that resemble those of planets orbiting 
slow-rotating stars so closely that the two cannot be dis- 
tinguished based on lightcurves alone. However, fitting 
the lightcurvcs of the planets around fast-rotators with a 
spherical-star model yields incorrect transit parameters. 
The significant systematic errors introduced into the deter- 
mination of the planet's radius are particularly distressing, 
especially at high impact parameters. This degeneracy can 
be resolved by independent determination of the stellar ro- 
tation period. 

Planets that arc not spin-orbit aligned can lead to spec- 
tacularly strange lightcurves that might not otherwise be 
immediately recognizable as planetary transits. Such plan- 
ets can yield asymmetric lightcurves and highly curved 
transit floors that rule out the spherical-star hypothesis. 
The precise shape of such transits reveals the stellar obliq- 
uity and the angle between the transit chord and the pro- 
jected stellar rotation axis. Together these two values de- 
termine the net spin-orbit alignment of the system. There 
is still a twofold degeneracy of the system, reflected around 
the projected stellar rotation axis. 

The changing temperatures across the face of the stel- 
lar disk lead to strong variability of transit lightcurves as a 
function of wavelength. Multi-wavelength transit photom- 
etry can then directly constrain the stellar temperatures, 
allowing for better models of the stars themselves. High- 
quality stellar models lead to more reliable measurements 
of the planet's radius. 

Perhaps ^ 5—10% of the Kepler target stars ought to be 
main sequence dwarfs of spectral type early F or earlier, 
and thus are probably rapid rotators. Determination of 
the spin-orbit alignment of planets orbiting these stars will 
provide a glimpse into the planet formation process around 
high-mass stars for the first time. 
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